Remarks on the renormalization of primordial cosmological perturbations 
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We briefly review the need to perform renormalization of inflationary perturbations to properly 
work out the physical power spectra. We also summarize the basis of (momentum-space) renormal- 
ization in curved spacetime and address several misconceptions found in recent literature on this 
subject. 
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I. INTRODUCTION 

The inflationary universe [l[ opens an exciting window 
to observationally test fundamental aspects of the theory 
of quantized fields in curved spacetime . In a curved 
background, the vacuum energy of quantum fields enters 
into the gravitational field equations. As is well known, 
infinities arise in the computation of vacuum energy, and 
other expectation values quadratic in fields, due to the 
ultraviolet behavior of the field theory. Therefore, the 
potential gravitational effects of the quantum vacuum 
must be handled with care. Methods have been devel- 
oped to define regularization and renormalization pro- 
cedures to physically account for the effects of vacuum 
energy in free and interacting quantum field theory in 
curved spacetime. 

In the cosmological scenario, it was shown in the early 
eighties that vacuum fluctuations can induce a primordial 
spectrum of density perturbations during an inflationary 
expansion. Remarkably, the calculated spectrum satis- 
factorily accounts for the origin of the cosmic inhomo- 
geneities that we observe in the present universe. Let <p 
represent a generic field describing (scalar or tensor) per- 
turbations during inflation. The quantum fluctuations of 
ip can be quantified by the mean square fluctuation in 
the vacuum state (see, for instance 
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where tp k (t) is defined by the expansion in Fourier modes 
of a free field operator 



<p(£,t) = / d 3 k(A, Vk (t) + A\<p* k (t)) ( 
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where and At are creation and annihilation opera- 
tors, such that Aj:\0) = 0. In the last equality of Eq. (JTJ 
we have defined the power spectrum A^,(fc,f), which is 
the quantity conventionally used in cosmology to quan- 
tify the variance (Q]) . For a single k, the power spectrum 
A^(fc,i) is well defined. However, the formal variance 
{tp 2 (x,t)) diverges in the ultraviolet. One could argue 
that this ultraviolet divergence concerns ultrashort wave- 
length modes, and it is not going to affect the predictions 
regarding the finite range of wavelengths that we actually 
observe in cosmic inhomogeneitics. But in view of the 
fact that the perturbations are a consequence of quan- 
tum field theory in curved spacetime, it was proposed in 
[§-[l2| that the physical power spectrum should be de- 
fined in terms of the renormaliz ed mean square 



(3) 



They found that, when methods of renormalization in 
curved space-time are applied, the power spectrum at 
wavelengths of observational interest may be affected in 
ways that can be tested by observations in the not-too- 
distant future. 

It has been argued recently that the renormalization 
procedure used in (cj-fl2| is not well-defined at the time 
scales at which cosmic inhomogenities are created dur- 
ing inflation, particularly at the time t k that each mode 
crosses the Hubble sphere [TH, Q] . In the present work, 
we briefly summarize the basis of renormalization (in- 
cluding its effects on the momentum-space power spec- 
trum) in an expanding universe and address the criti- 
cisms raised in |13|, Q] . We find that the arguments and 
conclusions proposed in [lU, [l4[ are in direct conflict with 
some of the basic principles of renormalization in curved 
space-time. The last paragraph in sections |TT] and IIIII 
summarize the arguments provided here in relation to 
the criticism of [lj, [l4[ • 
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II. RENORMALIZATION IN CURVED 
SPACE-TIME 

The renormalization of expectation values like those 
corresponding to the mean square fluctuation and the 
stress tensor of a quantized field in a curved space-time 
is more involved than in Minkowski space, even for the 
simplest case of a free field. This is rooted in the fact 
that in a general curved space-time there are insufficient 
isometrics to uniquely determine a global vacuum state. 
Moreover, the presence of space-time curvature yields 
new types of divergences not present in Minkowski space. 
This is true for free as well as interacting fields. 

It is instructive to briefly consider, in a general curved 
space-time, an interacting scalar field tp with a \tp in- 
teraction term in its Lagrangian. In flat space-time, one 
usually rcnormalizcs such an interacting theory by going 
to momentum-space and adding counter-terms to the La- 
grangian that will cancel the regularized UV divergences 
in the momentum space integrals corresponding to Feyn- 
man diagrams. The latter method is not directly applica- 
ble in a general curved space-time because of the absence 
of a global and generally-covariant momentum-space ex- 
pansion. However, this problem was overcome for a gen- 
eral smooth curved space-time in fl5j by making use of a 
Riemann normal coordinate system (RNC) with its origin 
at a given space-time point x' . Such a coordinate system, 
based on the system of geodesies leaving the point x' , ex- 
ists in any smooth space-time in a normal neighborhood 
of x' (i.e., one in which these outgoing geodesies do not 
intersect). In a RNC, the metric infinitcsimally close to 
x' is Minkowskian and has vanishing first derivatives with 
respect to the space and time coordinates of the RNC. In 
addition, the metric has a well defined expansion in pow- 
ers of the RNC coordinates y, which are defined for each 
space-time point x in the normal neighborhood of x' by 
means of the tangent vector at x' of the unique geodesic 
that connects x' to x and the invariant length of that 
geodesic. The coefficients in this expansion of the metric 
in powers of y are constants formed from contractions of 
powers of the Riemann tensor evaluated at the origin x' 
of the RNC system. 

In [l5|, Bunch and Parker (BP) defined a local 
momentum-space "Fourier" transform of the Feynman 
propagator G{x,x') based on an RNC with its origin at 
x' . Working in a general curved space-time, they used 
this local momentum-space method to evaluate the Feyn- 
man diagrams necessary to rcnormalize the interacting 
scalar field theory to second order in the interaction cou- 
pling constant A appearing in the Xtp 4 self-interaction. 
Dimensional regularization was used to replace the UV 
infinities that are present in 4-dimensions by covariant 
well-defined expressions. For the interacting theory, they 
introduced into the Lagrangian the minimal set of gen- 
erally covariant counter-terms necessary to absorb the 
regularized UV "infinities" into the values of the con- 
stant coefficients of these counter-terms. This process 
leaves one with a Lagrangian having terms of the same 



form as the original counter-terms, but with "renormal- 
ized' ' constant coefficients that are assumed to be finite, 
well-defined, and in principle measurable. 

The terms that involve the curvature at the space- 
time point x' were shown to drop out of the final result. 
This step is non-trivial and necessary for the renormal- 
ized interacting theory to be covariant in a general curved 
space-time. This calculation showed that interacting \tp A 
theory is renormalizable in a general curved space-time 
to second order in A. The reason for using the minimal 
set of counter-terms is to alter the original Lagrangian as 
little as possible, thus avoiding the arbitrary introduction 
of interactions that are not necessary for renormalization. 

They also calculated the leading terms of their lo- 
cal momentum-space expansion of the propagator that 
remain when A = 0, (i.e., for the non- interacting, 
free field). For example, for a minimally coupled free 
field, tp, the local momentum-space Fourier transform of 
\g(x)\~ 1 ^ 4 G(x, x')\g(x')\^ 1 ' 4 = G(x,x') has an asymp- 
totic expansion for large k 2 that is given by 

> (4) 

where k is the 4-momentum, m is the mass of the field, 
and the ellipsis (• • • ) includes terms that go as the third 
and fourth powers of l/(fc 2 + m 2 ). The coefficients of 
those terms are formed from contractions of Riemann 
tensors and their derivatives, evaluated at point x' , and 
arc given explicitly in Eq. (2.21) of [15j . The Feynman 
propagator in the normal neighborhood of x' is defined 
as in Minkowski space by replacing m 2 by m 2 — ie, with 
e small and positive, and carrying out the fco integration 
along the real axis in the complex plane. 

By carrying out the inverse Fourier transformation 
back to curved space-time, they were able to recast this 
series into the form of the proper-time or heat kernel ex- 
pansion of G(x,x'). They showed that the coefficients 
that they had calculated (which contain all terms up 
to fourth order in derivatives of the metric) in the se- 
ries in Eq.(j4]), when evaluated in the coincidence limit 
(x — > x'), are identical to the corresponding coefficients 
of the proper-time series that had been calculated [l6| 
for the free field in the coincidence limit. 

As is well known, the expectation values (T^{x')) of 
the energy-momentum tensor of the free scalar field in 
physical states can be found from the exact expression for 
the Green function G(x,x') by applying to G(x,x') the 
appropriate second-order differential operator that pro- 
duces the expression for the operator T^^a;') in the coin- 
cidence limit as x — > x' and taking the expectation value 
in the physical state of interest. From the asymptotic se- 
ries for G{k; x') discussed above, it is clear that the terms 
in that series up to those that go as l/(k 2 + m 2 ) 3 , will 
produce divergent contributions to (T^^x')} when the 
Riemann normal coordinate y is taken to zero and the 
dimension n is taken to the value 4 in the n-dimcnsional 
Fourier transform. It can be seen that, when the two spa- 
tial derivatives in T^ v arc taken into account, the first 



3 



term in ((4J giv cs a contribution to (T M1/ (a;')) that has 
a quartic UV divergence that does not depend on the 
curvature at x', the second term has a quadratic UV 
divergence that depends on the Ricci scalar curvature 
R(x'), and the subsequent terms have logarithmic UV 
divergences that depend in a more complicated way on 
contractions and spatial derivatives of the Riemann ten- 
sor, evaluated at x' . 

It is also clear that these UV divergent contributions 
to (Tfj, v (x')) have the same state- independent expression 
in terms of the curvature tensor and its derivatives at x'. 
The leading (quartic) divergence reduces to the vacuum 
energy in flat space-time and is customarily subtracted 
from the physical or renormalized value of (T^(x')) in a 
general curved space-time. The other state-independent 
divergences coming from the leading terms in the asymp- 
totic series in ((4]), are also subtracted. Alternatively, one 
can rewrite the full Lagrangian (including the gravita- 
tional terms) to include counter-terms having the same 
form as the expressions involving the Riemann tensor and 
its derivatives in these UV divergences. Then one can 
absorb the regularized expressions for these UV infinities 
into the constants in front of similar terms added to the 
original Lagrangian (similar to the procedure described 
above for the i nteracting field). Thus, the renormaliza- 
tion of the free field requires covariant terms involving 
the curvature tensor to be added to the Lagrangian with 
coupling constants that are in principle measurable. 

Whether one thinks simply of subtracting these state- 
independent UV-infinite terms from the formal expres- 
sion for T^y or thinks of those terms as renormalizing 
the coupling constants of terms present in the full La- 
grangian, the result is the same, namely, a quantity that 
has no UV infinities when its expectation value is eval- 
uated for any physical state in 4-dimensions. Its expec- 
tation value is what we will refer to as the renormalized 
expectation value of the energy-momentum tensor, de- 

noted by (T^(x'))re n -IH 

Now let us consider the renormalization of the vari- 
ance ((p 2 (x, t)) within the local momentum space method 
we are considering. The expression for the renormalized 
variance (<p 2 (x', t')) Ten in any given physical state is ob- 
tained by subtracting the terms in the asymptotic se- 
ries that would give UV infinities from the solution for 
G(fc;x') in the normal neighborhood of x' . The variance 
involves no derivatives of y(x), so only the first two terms 
in Eq. Q are subtracted from the local Fourier transform 
of G(fc;x'). Then the inverse of the local momentum 
space Fourier transform can be performed, giving a well- 
defined operator that we will denote by G(x; x') ron - The 
coincidence limit (x — > x') of the expectation value of this 
quantity, then gives the renormalized value of the vari- 
ance at x' in any given physical state. Thus, we obtain 
the result, 

(<£ 2 (x')>ren= lim (G(x; x')rc„) (5) 
x— >x f 

As explained previously, we are subtracting only the min- 
imum number of terms needed to give an expression hav- 



ing no UV divergences. Even though we found it nec- 
essary to subtract four terms of the asymptotic series to 
obtain (T^^x')),^, since the renormalization of the vari- 
ance involves only the first two counter-terms, there is no 
reason to subtract additional terms in the series. In the 
gravitational part of the Lagrangian of the renormalized 
theory the additional covariant terms needed to renor- 
malize the energy-momentum tensor will still appear, but 
it would be uneconomical and would lead to unnecessary 
complications to subtract the corresponding momentum 
space terms from G(k;x') in calculating the renormalized 
variance because those terms are not necessary to regu- 
larize UV divergences in the variance. One does not wish 
to introduce physical effects into the theory that are not 
required by the actual renormalization process. Henc e, 
one uses the principle of minimal subtraction. 

In the above discussion of renormalization in a general 
curved space-time, the asymptotic expansion of G(k; x') 
in (0| was used to identify the set of UV divergent terms 
in the local momentum space expansion in the normal 
neighborhood of x' . The asymptotic expansion was not 
used to approximate a solution of the differential equa- 
tion for G(fc; x') that would be the local momentum space 
Fourier transform of some particular global Green func- 
tion G(x,x') for values of x in the normal neighborhood 
of x'. Therefore, employing the asymptotic series up to 
fourth (or higher) order to rcnormalize the variance, ex- 
pecting that it will give a better approximation to the 
renormalized value, as unfortunately done in some liter- 
ature, is simply not correct. Similarly, it is incorrect to 
argue that the renormalization procedure fails for those 
values of k and R for which the second term in the ex- 
pansion Q is not small compared to the first term, as 
claimed in [lH E3] • This argument would spoil the full 
renormalization program. 



III. ADIABATIC REGULARIZATION IN AN 
RW UNIVERSE 

A spatially flat RW universe, for which ds 2 — dt 2 — 
a(t) 2 (dx 2 + dy 2 + dz 2 ), is often used as the background 
space-time in discussions of inflation. In such a universe, 
the 3-dimensional hypersufaces of constant t are homoge- 
neous and the relevant quantized perturbation field can 
be expanded as in Eq.Q. The ordinary differential equa- 
tion satisfied by the mode functions ifik(t) with suitably 
chosen boundary conditions can be solved exactly for cer- 
tain forms of the scale factor a(t). In such cases, one 
can also obtain the exact expression for the Green func- 
tion G(x,x') on any given spatial hypersurface. [2lj The 
method of adiabatic regularization in the RW universe 
0, [I?} starts with the formal expression for (T lw (x')), 
which is a function of t' in the RW coordinate system. As 
described above in connection with the local momentum 
space method, one can make a large momentum asymp- 
totic expansion of the 3-dimensional Fourier transform 
of this quantity and identify the leading terms in that 
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expansion that would give UV divergences in the inte- 
gration over k. These terms are the same for any phys- 
ical state. Now, instead of using dimensional regular- 
ization or some other regularization method to find the 
covariant form of the counter-terms that would be added 
to the Lagrangian in the course of renormalization, in 
adiabatic regularization one goes directly to the step of 
subtracting the relevant leading terms in this momen- 
tum space asymptotic expansion from the formal expres- 
sion for the spatial 3-dimensional Fourier transform of 
(T^x 1 )). This step insures that the integrand of this 
3-dimcnsional Fourier transform will not give any UV di- 
vergences when integrated over 3-momentum. The renor- 
malized expression for (T' All/ (x')) ren is defined by this mo- 
mentum integral. 

It is worth remarking that just as explained earlier in 
connection with the X(f> 4 interaction, the minimal num- 
ber of terms in the asymptotic series are subtracted in 
the process of adiabatic subtraction and renormalization, 
so as to change the form of the original unrenormalized 
Lagrangian as little as possible in constructing the renor- 
malized Lagrangian (including the gravitational part of 
the Lagrangian). 

Note that after making the subtractions in the above 
method of regularization the integrand of the momentum 
integral is already regularized (i.e., the momentum space 
integrand no longer has any UV infinities when it is con- 
sidered as a whole). If one wants to study the form of 
the counter-terms that would appear in the renormalized 
Lagrangian, then one could dimcnsionally regularize the 
individual subtraction terms and see what covariant cur- 
vature terms they would involve. However, in getting the 
renormalized result in the RW universe, there is no fur- 
ther regularization needed, apart from making the mini- 
mum number of subtractions we have described. One of 
the nice points of this process of adiabatic regularization 
is that it directly displays the spectral properties of the 
renormalized physical quantity, such as (T M y(x')) r en- 

If one is going to subtract the terms in the asymptotic 
series for all momenta, even small ones, then one has to 
resolve some ambiguities that are not determined sim- 
ply by the large k form of those terms. These possible 
ambiguities in the asymptotic form of the leading terms 
(such as whether or not to include the mass m 2 along 
with k 2 in the asymptotic series for large k), are resolved 
by requiring that during any intervals of time for which 
a(t) approaches a constant, the mode functions should 
approach linear combinations of the positive and nega- 
tive frequency forms that they would have in Minkowski 
space. In addition, one can require that in the limit of 
infinitely slowly changing a(t) the mode functions ipk{t) 
should take the form of the positive frequency Liouville 
(or WKB) adiabatic approximation to the mode function 
solution. This enforces the physically reasonable require- 
ment that the expectation value of the particle number 
operator is an adiabati c invariant; remaining unchanged 
in the limit of an infinitely slow and smooth expansion 
of the universe, regardless of the overall change in a(t) 



during the expansion. 

As outlined above for the energy-momentum tensor, 
the method of adiabatic regularization of the variance 
{tf 2 (x')) starts with the formal expression for (ip 2 (x')), 
which is a function of t' in the RW coordinate system. As 
described above in connection with the local momentum 
space method, one can make a large momentum asymp- 
totic expansion of the 3-dimensional Fourier transform 
of this quantity and identify the leading terms in that 
expansion that would give UV divergences in the inte- 
gration over k. These terms are the same for any physi- 
cal state. In adiabatic regularization one goes directly to 
the step of subtracting the relevant leading terms in this 
momentum space asymptotic expansion from the formal 
expression for the spatial 3-dimensional Fourier trans- 
form of (tp 2 (x')). This step insures that the integrand 
of this 3-dimcnsional Fourier transform will not give 
any UV divergences when in tegrated over 3-momcntum. 
The renormalized expression for (<p 2 (x')) rcn is defined by 
this momentum integral, regardless of the expansion rate 
of the universe (assuming the expansion is sufficiently 
smooth and that infrared divergences that may occur for 
zero mass are dealt with properly) . Therefore, it is again 
not correct to argue LT3l that adiabatic renormaliza- 
tion fails when the expansion rate H is larger than the 
physical momentum scale k/a(t), as already pointed out 
in the previous section. Additionally, in [14| a rather 
arbitrary redefinition of adiabatic subtraction is made, 
which docs not seem well founded. 



IV. RENORMALIZATION IN THE 
INFLATIONARY UNIVERSE 

In the physical situation in which the field tp(x, t) rep- 
resents scalar or (a polarization mode of) tensorial met- 
ric perturbation during inflation, ip(x,t) must be treated 
as a massless field. The term V" appearing in the wave 
equation of scalar perturbations, where V is the potential 
of a slow-roll inflationary model, should be considered as 
a second-order adiabatic term [l(| HH ■ In the massless 
limit, the subtraction terms defined by the BP renor- 
malization (l5j and those obtained with the adiabatic 
renormalization coincide, thus defining a unique expres- 
sion for the (renormalized) power spectrum [12| , namely, 
A 2 (/c,i) = 4-Kk 3 (\ip k (t)\ 2 -C k {t)). The subtraction terms 
Ck{t) can be obtained easily from the two first terms in 
the expansion Q (or equivalently by the adiabatic ex- 
pansion of the modes) . [22| 

For scalar perturbations in single-field, slow-roll infla- 
tion one finds (see [l2j for details) 
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2(2tt) 



a a o . 



5e) 



, (6) 



where Stf> is the inflaton perturbation (in the spatially flat 
gauge), a is the expansion factor, H = a/a , and e and 
r\ are the usual slow-roll parameters; while for tensorial 
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perturbations one gets 



1 



2(27r) 3 a 3 



a a . 



(7) 



where /i stands for one of the two independent fluctuating 
modes. The first term in square brackets in © and 
decays rapidly with time during inflation. However, the 
second term in square brackets varies very slowly dur- 
ing slow-roll inflation and, while maintaining the nearly 
scale-invariant behavior, it produces a non-negligible ef- 
fect in the power spectrum even at late times during in- 
flation, as explicitly worked out in (T^ |. 

Notice that the physical value of the correlation func- 
tion, call it (ip(x, t)ip{x' , i))ren, at two separated points 
x and x! , involves the subtraction terms in such a way 
as to insure continuity as x approaches x! . Rcnormal- 
ization can be naturally extended to higher-order point 
functions, as for instance the four-point function (odd 
point-functions are trivially zero) 



(ip{xi,t)ip(x 2 , t)<p(x 3 , t)ip(x4,,t)) Te 



= (ip(xi,t)ip(x 2 ,t)) lcn 

x (f(x3,t)(p(x 4 ,t)) Icn 

+ 2 perm. . (8) 



Therefore, it also maintains the Gaussianity of primor- 
dial scalar perturbations. However, the computations 
beyond leading order in cosmologial perturbation theory 
(producing non-Gaussianities) will be potentially affected 
by renormalization. Note that the simplest template 
to generate non-gaussianities [Tj| <p(x,t) = ipt{x,t) + 
fNL[<P 2 L (x,t) - (tp%(x,t)} Ien ], where tfL&t) denotes the 
linear Gaussian part of the perturbation, assumes implic- 
itly the necessity of renormalizing the variance. 
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There are some well known state-independent effects that 
can be calculated directly from the subtracted terms 
under discussion, as in the so-called conformal trace 
anomaly of the energy-momentum tensor. 
If one chooses in an RW universe any point x' on a 
given hypersurface at t' , then the RNC coordinate sys- 
tem with origin at x will have a normal neighborhood 
that includes the entire spatial hypersurface at t'. How- 
ever, these RNC coordinates are not the same as the RW 
coordinate system. 

In the massless limit the second subtraction term, see 
(j5J and (0, has an infrared logarithmic divergence. It 
can be naturally cured by introducing an infrared cutoff. 
We shall assume that this low-energy cutoff is far below 
the relevant cosmological scales. 



